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Abstract 

The aerodynamic loads on a Canard Rotor/Wing vehicle are investi- 
gated using inviscid numerical simulations in order to understand 
the flight characteristics of the vehicle during conversion from rotor 
craft to fixed-wing flight. A series of numerical simulations at seven 
azimuthal rotor indices are presented covering a quarter turn of the 
rotor. With symmetry arguments, these simulations produce 25 data 
points for a complete rotation. A Cartesian mesh approach is used 
to compute the flow field about a configuration with faired-over 
engine inlet and exhaust that matches the wind tunnel geometry. 
These simulations were performed using meshes with approximately 
9Million Cartesian cells. To better understand the aerodynamic 
effects of the rotor hub on the configuration, the same set of simula- 
tions were repeated for a hub-less geometry. Overall loads for both 
configurations are similar, but are due to somewhat different aero- 
dynamic mechanisms. 

1 Introduction 

The Canard Rotor/Wing (CRW) is a vertical takeoff and 
landing flight demonstration aircraft built via a joint Boeing/ 
DARPA effort. Figure 1 shows a 3-view of the configuration. 
The CRW is similar in concept to the X-wing prototype in 
that it uses a symmetric rotor with an elliptic airfoil section 
for takeoff, landing and hover, and that this rotor can be 
locked in a fixed position for use as a main-wing for high- 
speed forward flight. During hover, each rotor is cyclically 
pitched in order to generate lift; while in forward flight, both 
rotor blades are fixed in the wing position and collectively 
pitched to generate lift. 

During development, analysis and testing indicated that the 
X-wing design encountered high transient aerodynamic loads 
during the conversion from rotor craft to fixed-wing flight. 



Figure l: Three-view of the Canard Rotor/Wing flight demon- 
stration aircraft used in the numerical simulations. 


The CRW design addresses this issue by generating positive 
lift on both the canard and horizontal tail during conversion. 
By deploying the canard flap and the full-flying horizontal 
tail, enough lift can be produced to maintain altitude, the 
main rotor can be completely unloaded throughout the con- 
version maneuver. 

Despite the fact that the rotor is unloaded, it still decelerates 
through a non-uniform free stream, and therefore will still 
transmit unsteady loads to the airframe. The frequency of the 
loads will be the harmonics of the rotor speed, but then- 
amplitude is difficult to estimate. Prediction of the amplitude 
of these unsteady loads is essential in sizing the control sur- 
faces and actuators of the flight control system. 

In this work, we use numerical simulations to get quantitative 
estimates of the magnitude of these unsteady loads and an 
understanding of the aerodynamic mechanisms which pro- 
duce them. Work on the X-wing indicated that the most 
severe loadings occurred during the last few revolutions of 
the rotor when it turns relatively slowly. Under these condi- 
tions, the rotation time-scale is much smaller than that of the 
flow, making a sequence of steady measurements at fixed 
rotor azimuthal positions an efficient and viable method to 
measure and predict the loads. For this reason, the wind tun- 
nel experiments^ and numerical simuladons(in this paper) 
were carried out with the rotor fixed at a series of azimuthal 
positions around the rotor disk. 

The simulations are performed with the rotor positioned at 0°, 
15°, 30°, 45°, 60°, 75°, and 90°. Symmetry extends the data 
from these seven first quadrant rotor positions to fill 25 sta- 
tions for one 360° turn of the rotor. The flight conditions for 
conversion from rotor craft to fixed-wing flight are set at a 
Mach number of 0.1986(130Knots) and 0° angle-of-attack. 
Each of the computations is performed using the Cartesian 
mesh flow simulation package discussed in references [2-4], 
The surface preparation, triangulation, and Cartesian mesh 
generation capabilities in this package are used to generate 
the mesh system at each of the azimuthal positions. Two 
inviscid flow solvers are included in this package and both 
have been extensively validated in the open literature. 13 ' 

This paper first presents preliminary solutions obtained on 
the complete CRW with the hub included. The integrated 
forces and moments from these simulations are presented as 
a function of the rotor’s azimuthal position. The trends in 
these forces and moments are explained in terms of the 
changes in surface pressure. These simulations suggest that 
the hub has a significant effect on the flow and thus on the 
vehicle loads. In order to better understand the hemispherical 
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hub’s contribution to the changes in the flow and the loads, a 
full series of simulations are presented with the hub removed. 
Though this case is not realistic in that there is no hardware 
attaching the rotor to the fuselage, it helps us understand the 
aerodynamic mechanisms that contribute to the changes in 
vehicle loads. 

2 Methods 

The process begins with receiving geometry components as 
CAD solids. The components can be moved, and/or rotated 
with respect to each other in order to perform parametric 
studies. In the current case, the rotor is positioned at several 
azimuthal angles using this feature. Surface triangulations are 
generated for each of the CAD solids using a combination of 
CAD tools and auxiliary tools bundled with the Cartesian 
mesh flow simulation package. These component triangula- 
tions are then provided as input into the Cartesian mesh gen- 
eration system. 

After indexing the rotor to the desired azimuthal position, the 
component triangulations are intersected together to extract 
the wetted surface of the configuration. The intersection 
program in the Cartesian package[2] incorporates adaptive 
precision arithmetic, and includes an automatic tie-breaking 
routine that consistently resolves geometric degeneracies 
without user intervention. This tool computes the intersec- 
tions between closed body triangulations, and retriangulates 
the intersection neighborhood to join the pieces together into 
a single triangulation. 

The volume mesh generation algorithm takes the intersected 
surface triangulation as input. This program generates an 
unstructured Cartesian mesh through subdivision of hexahe- 
dral cells of an initially uniform coarse grid. It adapts to the 
contours of the geometry, thus creating a mesh where the size 
of the cubic cells is smallest in the neighborhood of the 
sharpest features of the geometry, while away from the 
geometry, the grid remains coarse. Specific regions of 
enhancement can also be defined in order to ensure that a 
given mesh size is maintained in a specified area. The pro- 
gram also computes the intersections of the geometry with 
the cells of the mesh which pierce the geometry in order to 
define polyhedral cut-cells. The faces of these polyhedral 
cells are then used to specify the tangency boundary condi- 
tion in the flow solver. Meshes generated by this algorithm at 
five different rotor indices are presented in Figure 2. The 
body and cuts from the mesh in the vicinity of the body are 
shown. Cuts through the mesh, behind and under the body, 
make up the background of the figures. One small figure 
shows the canard/flap region in greater detail. The approach 
in this volume mesh generator preserves the ability to refine 
the mesh to different degrees in different directions. This 
makes it possible to avoid generating an excessive number of 
Cartesian cells in three dimensions. 

A Cartesian flow solver that solves the Euler equations is 
used to compute the flow field. The Euler equations describe 
the unsteady flow of an inviscid perfect gas and may be 
expressed in integral form. 


a 

Here Q is the state vector of conserved quantities, and F is 
the tensor of flux density for an inviscid perfect gas. refers 
to an arbitrary control volume and 9Q is its closed boundary 
with an outward facing surface normal vector n . 

The flow solver is an upwind code that uses an explicit multi- 
stage Runge-Kutta time stepping algorithm to drive the solu- 
tion to a steady state. An optimally damped 2 nd order scheme 
using five Runge-Kutta stages was used to compute the 
present set of subsonic cases[6]. The flow solver also uses 
multi-grid convergence acceleration to damp the high fre- 
quency error modes. In the present computations, five levels 
of multi-grid are used. The coarse meshes for the multi-grid 
computations were generated by a utility in the Cartesian 
mesh flow simulation package. This utility is documented in 
reference [3]. Details of this flow solver scheme can be found 
in the references [3,4]. 

3 Simulations 

The solutions are presented below for all rotor settings at a 
Mach number of 0.1986 and 0° angle of attack. This corre- 
sponds to the conversion scenario where the aircraft is 
switching from hover to forward flight at 130 Knots. For all 
simulations, the boundaries of the mesh spanned 20 times the 
body length in all directions. 

The change in the value of lift is used to determine the con- 
vergence of the solution. All cases are run untill the value of 
lift converges to within 1%. Approximately two orders of 
magnitude drop in the residual on the finest mesh after a full 
multigrid startup is also required. 

3.1 Simulations with hub 

A simulation of the full geometry with the hub is conducted 
for seven rotor positions in the first quadrant. As discussed in 
the introduction, this results in 25 data points, due to symme- 
try, covering one 360 degree rotation. In this section, we dis- 
cuss the trends found during these simulations. 

Figure 3 shows the body of the CRW, colored by the coeffi- 
cient of pressure, for all computed rotor positions. In these 
surface pressure figures, high pressure corresponds to 
magenta, while suction corresponds to blue. The suction on 
the top surfaces of the canard and horizontal tail shows that 
these components are producing lift. This is seen clearly in 
Figure 5 which shows the componentwise breakdown of the 
lift. The rotor produces little or no lift in comparison with the 
canard and horizontal tail. In other words, the aircraft is fly- 
ing entirely on the lift produced by the canard/flap and the 
horizontal tail. The aircraft is no longer dependent on the 
rotor for lift. Figure 5 also shows that the lift produced by the 
fuselage, the canard, the horizontal tail and the vertical tail 
does not change very much as the rotor turns. The changes in 
total lift are almost entirely due to the changes in lift pro- 
duced by the rotor. Upon close inspection of Figure 3, we 
realize that most of the surface pressures do not change from 
rotor position to rotor position. However, the pressure on the 
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underside of the rotor seems to change with the changing While the changes in pressure are minimal over most of the 
rotor position as shown in Figure 4 on page 4. surface, in the region of the rotor hub fairing the surface pres- 

sure changes from rotor position to rotor position. A possible 
mechanism for changing the flow on the rotor is the down- 




grid is indicated. The close-up shown in the final frame is representative of the resolution in all meshes. 
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wash from the canard/flap. If the rotor sees the downwash 
from the canard/flap, it would produce negative lift simply 
because the effective angle of attack on the rotor would then 
be negative. Figure 6 shows streamlines in a vertical plane 
that cuts the outboard part of the canard/flap for the 90° rotor 
position case. This figure shows that while there is a down- 



Figure 3: Surface pressure on the CRW for seven rotor positioi 



Figure 4: Underview of the surface pressure 
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wash associated with the lift generated by the canard/flap, the 
rotor is not effected by this down wash. Thus, the changes in 
pressure on the underside of the canard must be due to other 
effects. 


CRW componentwise breaJoiown of C 1 vs. azimuth angle 



Having ruled out the downwash from the canard/flap, there 
are two other factors contributing to these changes in surface 
pressure under the rotor. One factor is the upwash from the 
front part of the fuselage. The upwash is a result of the fuse- 
lage pushing the flow upward on the front side of the vehicle. 
When die rotor is close to the fuselage, a large part of the 
rotor is exposed to this upwash generating forces and 
moments which diminish as the rotor moves away from the 
body. Another factor is the suction on the rotor hub fairing 
assembly. This strong suction is the result of the flow acceler- 
ating around the hub. From the surface pressure(Figure 4) we 
note that the peak suction follows the rotor’s azimuthal 
progress. In other words, the amount of suction seen by the 
underside of the rotor is larger as the rotor moves away from 
the body. 



Figure 6: Effect of the downwash from the canard/flap on the 
rotor 


The variations in the suction acting on the horizontal and ver- 
tical surfaces around the hub along with the effects of the 
upwash correlate with the changes in overall vehicle forces 
and moments. The trends of the flow can be best examined if 
we look at the variations in forces and moments with respect 
to the azimuthal position of the rotor and correlate them to 
the visible changes in pressure on the surface. The plots of 
the total forces and moments with respect to the azimuthal 
position of the rotor are presented in figures 7 through 12. 
These plots include the curves for the case without the hub, 
which are discussed in section 3.2 on page 7. Due to the pro- 
prietary nature of this data, the values can not be included on 
the y-axis of figures 7 through 12. 

A plot of the total lift with respect to the azimuthal angle in 
Figure 7 shows that more lift is being generated when the 
rotor is in the 0° position. Upwash from the fuselage is 
largely responsible for this higher lift. The upwash can 
clearly be seen in Figure 3 on the front of the ftiselage. In 
Figure 4, the underside of the rotor at 0° has higher pressure 
than when the rotor moves away from the body. This contrib- 
utes to more lift at 0® and the lift gradually reduces as the 
rotor advances since less and less of the rotor is in the 
upwash field. As the rotor moves away from the body, the 
suction on the sides of the hub and on the underside of the 
rotor becomes the dominant effect. As the rotor approaches 
the 90° position, the suction under the rotor gets larger. 
Figure 4 shows the surface pressure from the underside of the 
rotor so that the hub region is readily visible. We observe that 
the largest suction on the hub is created when the rotor is in 
the 90° position. This suction generates a negative lift on die 
rotor/hub assembly, which contributes to a lower total lift. 
Thus, the 90° case has lower lift than the 0° case. 

C { vs. Azimuth Angle; M=0. 1986, a=4), (i=0 



The variation in the drag force, shown in Figure 8, is much 
smaller and other effects such as the separated flow behind 
the hub and the faired-over engine exit may be the primary 
causes of drag variation. In order for the suction in the hub 
region to have a direct effect on the drag force, the suction 
must act on vertical surfaces normal to the flow direction in 
the vicinity of the hub. Such surfaces do not exist. Thus, the 
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movement of the suction in the hub region does not have a 
large effect on the drag. A possible source of drag is the inter- 
ference between the rotor and the fuselage. This interference 
is largest when the rotor is right on top of the fuselage (the 0° 
case). The interference drag, therefore is largest at the 0° 
rotor position and reduces as the rotor moves away from the 
body. 

C 4 vs. Azimuth Angle; M -0.1986, Ct=0, ^=0 



The side force, shown in Figure 9, can also be explained by 
the tendency of the suction to follow the rotor. In this case, 
when the rotor is in the 0° or the 90° positions, the flow is 
symmetric and the side force is zero. However, in all other 
cases, the side force is non zero due to the fact that the suc- 
tion is higher on one side of the body than the other. This dif- 
ferential suction acts on the vertical surfaces close to the hub, 
thus generating a side force. The sign of the side force can 
also be explained by the suction from the hub. If we look at 
the 45° case in Figure 4, we see a larger amount of suction on 
the back rotor as compared to the suction under the front 
rotor. The result is a positive side force. This is confirmed in 
Figure 9 where the side force is indeed positive at the 45° 
rotor position. The 45° case should be the same as the 135° 
case, but with the opposite effects since the rotor is on the 
other side of the body. Thus, the value of the side force for 
the 135° case should be negative of the 45° case. The logic 
carries over to the rest of the cases between 90° and 180° and 
again between 270° and 360°. The values between 180° and 
270° are duplicates of the values between 0° and 90°. 

The variation of the moments with respect to the azimuthal 
angle of the rotor can also be explained in terms of the suc- 
tion created by the hub. The suction on the hub follows the 
rotor; thus at a rotor position of 45°, the suction is greater 
close to where the rotor is. However, due to the oncoming 
flow, the suction in the front is relieved and so closer to the 
front facing rotor, there is less suction than the back side. 
This contributes to a larger downward force on the rearward 
swept rotor which in turn contributes to a more positive 
pitching moment. In addition, the upward sloped fuselage 
induces an upwash on the rotor as evidenced by the higher 
pressure on the bottom of the front rotor close to the fuselage. 
This further contributes a positive pitching moment. Because 


vs. Azimuth Angle; M=0. 1986, ct=0, j}=0 



the effect is seen by more of the rotor when the rotor is closer 
to the body, the moment is higher close to the 0° position than 
when most of the rotor is away from the body at the 90°posi- 
tion. The total variation in pitching moment can be controlled 
with approximately 2° of elevon deflection. 


Cpuck v * Azimuth Angle: M -0.1986, a=43, fH) 



The rolling moment, though much smaller in magnitude than 
the pitching moment, also depends on the hub suction behav- 
ior. Once again, we compare the 0°, 45° and 90° cases. At 0 
and 90 degrees, the rolling moment is zero due to the symme- 
try of the flow. However, at 45°, the suction has clearly 
caused the right rotor to produce less lift than the left rotor. 
For this reason, there is a positive rolling moment produced. 
The upwash also contributes to a higher lift on the left rotor. 
This is ascertained from Figure 4 as there is a high pressure 
region under the left rotor on the side closer to the ftiselage. 
Similarly, a positive rolling moment is produced for the 15°, 
30°, 60°, and 75° cases. Like the side force, the symmetry of 
the flow dictates the sign of the moment at other rotor posi- 
tions. The total control authority of the elevon is also marked. 

A small yawing moment is produced due to the suction force 
acting differentially on the vertical surfaces close to the hub 
and its sign can also be explained in terms of the hub suction. 


6 - 







American Institute of Aeronautics and Astronautics, AIAA Paper 2001-0997 


vs. Azimuth Angle; M =0.1986, O»0, fM) 



The yawing moment can be countered with a 3° deflection of 
the rudder. 


vs. Azimuth Angle; Si=0. 1986, a=0, (M) 



3.2 Simulations without hub 

To isolate the aerodynamic effects of the hub, a simulation at 
the same rotor positions as the previous case is conducted 
without the hub. This is not a realistic configuration since the 
mast and control arms are not included in the geometry. 
However, it represents the extreme case and helps understand 
the effect of the hub on the force and moment trends. 

The results are presented in Figure 1 3 with surface pressure 
and in figures 7 through 12 with the line plots of the forces 
and moments with respect to the azimuthal angle. Here also, 
we look to the surface pressure to explain the trends in the 
forces and moments. Two reasons to expect differences 
between the flows over the geometry with the hub and the 
one without are evident. First, taking the hub away removes 
the blockage and thus the mechanism for the flow to acceler- 
ate around the blockage under the rotor. Second, a channel 
with top and bottom walls is created with open sides through 
which flow will accelerate due to the venturi effect. The ven- 
turi effect creates a suction that is less severe but more con- 
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centrated under the center of the rotor than the case with the 
hub. The trends in forces and moments are very similar to the 
case with the hub as seen in figures 7 through 12. We use sur- 
face pressure shown in Figure 13 to explain these trends. 

The trend in lift is very similar to the case with the hub 
because just as there was a larger region of low pressure 
under the rotor for the 90° case, so it is true for the case with- 
out the hub. The surface pressure seen from the underside 
shows clearly that while the suction due to the accelerated 
flow is of much smaller magnitude in the 0° case, it is the 
dominant feature for the 90° case. So the rotor generates 
more negative lift when the rotor tips are away from the body. 
In addition, the upwash from the fuselage is also present and 
adds to the higher lift for the cases when the rotor is in the 
upwash field. We can also explain why the lift is greater for 
the case with the hub than the case without the hub. The hub 
has a lifting component in the front and in the rear where the 
pressure is high due to the fact that the flow is being blocked 
by the hub. This high pressure contributes to lowering the 
amount of negative lift contributed by the rotor/hub assem- 
bly. The case without the hub has a suction under the rotor, 
but no high pressure area to counter it. So the overall lift of 
the case with the hub is larger than the one without the hub. 

The changes in drag follow the same trends as the case with 
the hub for the same reasons. This time, however, the hub is 
not present which results in less blockage and thus in a lower 
overall drag. 

The side force once again has the same trend, but it is lower 
than the case with the hub. This is because the flow which 
spilled over onto the sides of the fuselage from the hub creat- 
ing a suction region on the vertical walls of the fuselage is no 
longer present. The non-symmetric effect of the suction 
under the rotor on the flow on the vertical walls of the fuse- 
lage is not as pronounced, resulting in lower side force. 

The pitching moment trend is also similar to the case with the 
hub. Once again we see that the venturi-suction creates a 
larger suction region on the back rotor and the same argu- 
ments we made in the hub case work here. The high pressure 
due to the upwash from the fuselage on the part of the left 
rotor close to the fuselage is also contributing positive pitch- 
ing moment. 

The rolling moment is also the victim of the differential suc- 
tion under the rotor just as the case with the hub is. The 
upwash field responsible for the higher lift on the left rotor 
also contributes to the rolling moment as it did for the case 
with the hub. The actual magnitude is one order smaller than 
the pitching moment and the difference between the 
“hubbed” and “hubless” cases is minimal. The upwash is the 
same in both cases, but the differential suction due to the hub 
has now been replaced by the differential suction due to the 
venturi effect. 

The yawing moment can also be explained in the same man- 
ner and it is also approximately the same as the hub case. 


4 Concluding Remarks 

The vehicle loads on the CRW are evaluated during the con- 
version from rotor craft to fixed-wing flight using computa- 
tional tools for inviscid analysis with a Cartesian mesh. The 
resulting trends in forces and moments are explained in terms 
of the forces imparted on the body due to the aerodynamic 
pressure. 

It is confirmed that the rotor is indeed unloaded or close to 
unloaded during conversion. It is also confirmed that the load 
of the aircraft can be carried by the canard and the horizontal 
tail. Plots of the variation of the forces and moments with 
respect to the rotor position are generated and explained in 
terms of the varying suction in the rotor hub fairing region for 
the case with the hub and in terms of the varying venturi-suc- 
tion in the channel between the rotor and the fuselage for the 
case without the hub. The upwash from the upward slanting 
fuselage is also responsible in part for the changes in loads 
with respect to the indexing of the rotor. 

The investigation of the results suggests that the case of a 
hemispherical hub imparts similar loads to the vehicle as the 
case where the hub has been removed. Though the mecha- 
nism for the flow acceleration under the rotor is different in 
each case, the effective changes in forces and moments are 
quite similar. The contribution of the upwash field from the 
front part of the fuselage to the loads on the rotor remains 
unchanged regardless of the presence of the hub. 

The simulations and the investigation of the cause of the vari- 
ation in the forces and moments suggest that the forces and 
moments imparted to the geometry may be substantially dif- 
ferent when the mast and other hardware are present 
between the rotor hub and the upper fuselage. This will 
relieve the suction under the rotor while introducing some 
blockage in the flow. No changes in the geometry of the hub 
area can relive the upwash from the forward fuselage. Thus, 
the effect of the upwash is the limiting case for the load vari- 
ations experienced by this vehicle. Further study is required 
to ascertain the changes in loads due to this kind of modifica- 
tion to the geometry. 

A non-hemispherical, “more aerodynamic” hub geometry, 
rather than the complete removal of the hub may also reduce 
the unsteady loads imparted to the vehicle. However, we also 
leave this for future investigation. 
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